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ABSTRACT 

The Galactic disk retains a vast amount of information about how it came to be, and how it evolved 
over cosmic time. However, we know very little about the secular processes associated with disk 
evolution. One major uncertainty is the extent to which stars migrate radially through the disk, 
thereby washing out signatures of their past (e.g. birth sites). Recent theoretical work finds that such 
"blurring" of the disk can be important if spiral arms are transient phenomena. Here we describe an 
experiment to determine the importance of diffusion from the Solar circle with cosmic time. Consider 
a star cluster that has been placed into a differentially rotating, stellar fluid. We show that all 
clusters up to ~ 10 4 M© in mass, and a significant fraction of those up to ~ 10 5 Mq, are expected 
to be chemically homogeneous, and that clusters of this size can be assigned a unique "chemical 
tag" by measuring the abundances of < 10 independent element groups, with better age and orbit 
determinations allowing fewer abundance measurements. The star cluster therefore acts like a "tracer 
dye" , and the present-day distribution of its stars provides a strong constraint on the rate of radial 
diffusion or migration in the Galactic disk. A star cluster of particular interest for this application is 
the "Solar family" — the stars that were born with the Sun. If we were able to identify a significant 
fraction of these, we could determine whether the Sun was born at its present radius or much further 
in. We show all-sky projections for the Solar Family under different dynamical scenarios and identify 
the most advantageous fields on the sky for the experiment. Sellwood & Binney have argued for strong 
radial transport driven by transient spiral perturbations: in principle, we could measure the strength 
of this migration directly. In searching for the Solar family, we would also identify many thousands of 
other chemically homogeneous groups, providing a wealth of additional information. We discuss this 
prospect in the context of the upcoming HERMES high-resolution million-star survey. 
Subject headings: Galaxy - stellar populations - open clusters elemental abundances 



1. INTRODUCTION 

An axiom of chemical evolution models is that stars 
form in situ within the disk, or at least they have not 
migrated far from their formation radius. But it has long 
been recognized that the metal abundance scatter at all 
radii is very large and comparable to the amplitude of the 
overall declining slope in [Fe/H]. This led to the idea that 
either the ISM was much less homogeneous in the past or 
that stars have migrated over large radial distances (Wie- 
len, Fuchs & Dettbarn 1996). The allowed distribution 
of orbit ellipticities accounts for only half of the scatter 
(e.g. Nordstrom et al 2004). Radial migration requires 
a strong non-axisymmetric impulse, i.e. perturbations 
from massive clouds or spiral arms. Cloud deflections are 
closely related to spiral arm impulses because the orbit- 
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ing cloud generates a spiral response ( "swing amplifier" ) 
within the disk (Julian & Toomre 1966). In a seminal 
paper, Sellwood & Binney (2002; hereafter SB02) went 
on to show that transient spiral perturbations can lead 
to strong radial "churning" in disk galaxies (i.e. mixing 
without heating), and this can reasonably explain the 
observed abundance scatter (Schonrich & Binney 2009). 

In a frame that rotates with the pattern speed of the 
spiral arm, stars close to the co-rotation radius orbit a 
point moving at the pattern speed. Inside co-rotation, 
stars overtake the pattern and fall into the spiral po- 
tential; this boosts their angular momentum (measured 
about the rotation axis of the disk) which moves them to 
a larger radius. The ejected star now lags behind the pat- 
tern speed and falls backwards into the wave approach- 
ing from behind. This lowers its angular momentum and 
the star returns from whence it came which results in no 
overall migration for an individual star. The key insight 
is to consider what happens if spiral arms are transient 
(e.g. Barbanis & Woltjer 1967). This results in a strong 
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stellar migration with individual stars moving inwards 
and other stars moving outwards - the Sellwood-Binncy 
mechanism. Stars near resonance can lose or gain a sub- 
stantial fraction of their angular momentum in a single 
encounter. This process has now been seen in N-body 
simulations (gas+stars) of evolving disks that are kept 
sufficiently cool from halo gas accretion (Roskar et al 
2008; Sanchez-Blazquez et al 2009), and the radial churn- 
ing appears to be strong. 

In support of radial migration, the Sun appears to be 
+0.17 dex too metal rich for its environs (Edvardsson et 
al 1993) which led Wielen et al (1996) to suggest that 
the Sun has moved outwards by about 2 kpc over the 
course of its lifetime. SB02 note that the Sun's highly 
circularised orbit is consistent with having migrated to 
its present position through churning. All stellar popula- 
tions, from the oldest (e.g. planetary nebulae; Maciel et 
al 2006) to the youngest (e.g. Cepheids; Andrievsky et 
al 2004), exhibit a marked abundance gradient and show 
evidence of having flattened with cosmic time. This pro- 
vides some justification for associating a metal-rich Sun 
(cf. Asplund et al 2005; Haywood 2008) with a smaller 
birth radius. Interestingly, Roskar et al (2008) find that 
metal-rich stars, like the Sun, could have originated al- 
most anywhere. They find that half of the stars within 
the Solar torus today have come from large radial dis- 
tances (>2 kpc). 

So where are the stars that were born with the Sun? 
If we were able to identify a useful fraction of the Solar 
familyj^] we would be able to answer this question. We 
can think of the Solar family as a tracer dye that has 
been dropped into a differentially rotating, stellar fluid. 
If the Sun was born at its present radius, in the absence 
of strong radial scattering, we would assume that the So- 
lar family is confined to the Solar torus. In this picture, 
the width of the annulus is likely to reflect the epicyclic 
radial amplitude (~ 1 kpc) of stars with the age of the 
Sun (§5.3). If the stars are found to lie along a thin- 
ner annulus, this would indicate possibly that the orbits 
have been confined by strong resonances from the cen- 
tral bar, say. If most of the Solar family was found to 
lie inside of the Solar circle, we would infer that it has 
experienced a steady outward radial migration. Grenon 
(1999) considered the possibility that stars are formed 
near the Galactic centre and then diffuse outwards. If 
we were unable to find any siblings, this may present the 
biggest surprise of all (see §2). 

The important issues of birth place and radial migra- 
tion are difficult to unravel. Just how secular migration 
evolves depends in large part on the Galaxy's accretion 
history which sustains and triggers the onset of spiral dis- 
turbances (Sellwood & Carlberg 1984). The stars that 
are most affected by angular momentum transfer lie close 
to the corotation radius, but the stochastic nature of the 
spiral disturbances give rise to a range of pattern speeds 
over the same time interval (Sellwood 2008). In partic- 
ular, in the SB02 simulation, the power spectrum of the 
m = 2 modes reveals a wide range of pattern speeds 

3 The Solar family was first discussed by Bland-Hawthorn &c 
Freeman (2004): these siblings should not be confused with "Solar 
twins" (e.g. Hardorp 1978; Dorren &; Guinan 1994; Cayrel de 
Strobel 1996), i.e. stars that are essentially identical in character 
to the Sun but at different stages of their evolution. By definition, 
these cannot be members of the Solar family. 



(their Fig. 10), with corotation radii in the range 4 kpc 
to 20 kpc. 

If the Sellwood-Binney mechanisrrF] is as efficient as 
suggested by Roskar's models, we are forced to conclude 
that the conventional approach to chemical evolution 
modelling (cf. Chiappini et al 2001) must be overhauled 
(Schonrich & Binney 2009). Furthermore, we must em- 
brace a broader range of chemical elements in surveys in- 
volving millions of stars (see Bland-Hawthorn & Freeman 
2004). Most stars are born in a single burst within com- 
pact clusters and stellar fragments. We can use chemical 
tagging to ascertain whether these substructures are ho- 
mogeneous systems or not. De Silva et al (2006, 2007a, b) 
have recently demonstrated that open clusters are chem- 
ically homogeneous, and even moving groups can show 
the same signature (see also Castro et al 1999; Sestito 
et al 2007; Randich et al 2006; Shen et al 2005; Chou 
et al 2010). Apart from a few light elements, globular 
clusters are also highly homogeneous (Gratton, Sneden 
&; Carretta 2004). The fact that old clusters are chemi- 
cally homogeneous argues strongly against internal stel- 
lar processes or external ISM accretion as influencing the 
chemical signatures. We can therefore reconstruct many 
stellar clusters through their unique chemical signatures. 

It is now possible to obtain echelle spectroscopy (R > 
30, 000) for millions of stars in campaigns of a few 
years. The technique of "chemical tagging" (Freeman 
& Bland-Hawthorn 2002) is being addressed by ongoing 
and planned stellar surveys. These include the HERMES 
survey on the Anglo- Australian Telescope (Freeman & 
Bland-Hawthorn 2008; Freeman, Bland-Hawthorn & 
Barden 2010), and the APOGEE survey at Apache Point 
Observatory (Allende Prieto et al 2008). 

The structure of our paper is as follows. In §2, we 
estimate the birthrate of clusters within the Solar torus 
as a function of mass. In §3, we examine their likely 
chemical uniformity as a function of the progenitor cloud 
mass. In §4, we derive more stringent limits on the size 
of the Solar family at the time of its birth. In §5, we 
simulate four possible distributions for the Solar family 
and determine their projection on the sky. In §6, we 
describe an experiment aimed at finding members of the 
Solar family before a discussion of the implications in §7. 

2. FORMATION OF STAR CLUSTERS 

If we consider our volume of interest to be the Solar 
torus, defined as the annulus 2 kpc thick in r adiu s cen- 
tered on the Solar circle at R « 8 kpc (see 
the cluster birth rate within this region is 



5.1 1, then 
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dtd\n{M.,/M Q ) 



3 x 10 3 M 
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where Mq — M/(M© yr *) and M is the total star 
formation rate in the Galaxy and the upper and lower 
cluster mass limits are (M*, m i n , M» >max ) = (5OAf ,2 x 
10 5 M Q ). Thus, for example, in the last Gyr we ex- 
pect ~ 10 4 clusters whose masses are within a factor 
of three of 10 4 Mq to have formed within the Solar torus 

4 These models treat only stellar migration. It is interesting to 
consider whether gas experiences the same churning or whether 
the shor t lifetime of giant molecular clouds (~ 30 Myr; |Fukui et] 
|al.|2009[ l suppresses it. 
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if the star formation rate over this time has been close 
to its present value of ~ 3 Mq yr -1 (McKee & Williams 
19971. This result is quite insensitive to the assumed 
limits on the cluster mass function or the extent of the 
star- forming disk. Changing M 4 m i„ or M +max by a fac- 
tor of 10 alters the estimated cluster birth rate by only 
~ 20%. 

All but a small fraction of clusters formed are unbound, 
and the stars born within them drift off to become part 
of the field star population. By comparing catalo gs of 
open and embedded clusters, Lada fe Ladaj ( |2003[ ) esti- 
mate that less than 10% of clusters survive more than 10 
Myr, and that only 4 - 7% survive for 100 MyiQ More- 
over, analysis of extragalactic cluster populations (albeit 
with masses significantly larger than the typical Milky 
Way cluster) suggests that the survival fraction is nearly 
independent of mass until two-body relaxation becomes 
the dominant destructi on process at times ^ > 100 Myr 
after cluster formatio n (|Fall et al.|2005[ |2009[ [FaiT||2006 



Whitmore et al. 20071. During this first 100 Myr, roughly 
80 — 90% of clusters are disrupted independent of mass, 
while two-body relaxation at longer times preferentially 
destroys low mass clusters. Since we do not strive for bet- 
ter than 10% accuracy in our prediction of the disrupted 
cluster birth rate, we will not introduce the added com- 
plication of a mass- or time-dependent disruption rate. 
Instead we simply assume that a fixed 95% of all clusters 
that are born are disrupted, independent of mass and 
time. This will slightly overestimate the number of dis- 
rupted clusters that formed less than 100 Myr, but for 
a relatively quiescent star formation history such as the 
Galaxy's, these clusters constitute a negligible fraction 
of the total population. Thus, if we target stars older 
than 100 Myr, the birth rate of now-disrupted clusters is 
given by equation [I] to within 10%. 

3. CLUSTER CHEMISTRY 
3.1. Introduction 

Since most stars are not born in clusters that remain 
bound, and thus do not end up in open clusters more than 
10 Myr later, the clusters observed by De Silva and col- 
laborators presumably represent a biased sample. They 
are either the innermost remnants of clusters that lost 
most of t heir members and pro duced only a small bound 
core (e.g. Kroupa et al. |[2001 1 , or they are clusters that 
remained bound because they had star formation effi- 
ciencies significantly higher than the typical 10 — 30% 
(Lada & Lada |2003| ) . We would therefore like to define 
some rough theoretical expectations about the chemical 
homogeneity of the disrupted cluster population. 

Mixing of heavy elements in protoclusters is governed 
by several time scales. The dominant mode of chem- 
ical mixing in molecu lar clouds is turbulent diffusion 



(Murray & Lin 1990). The t ime scale associated with 
H 2 /K flXie et al. 



1995), where 



this process is idiff ~ 

H = (ifr — H^)^ 1 is a composite between the scale 

heights of the H 2 molecules, H^ 2 = ™H 2 /|VnH 2 |, and 
of the species i whose diffusion time we are computing, 
Hi = rij/|Vnj|. Here nn 2 and are the number densi- 
ties of the two species. The quantity K is the diffusion 

5 We point out, however, that the distribution must have a long 
tail because a few open clusters have ages 2> 1 Gyr (Friel 1995). 



coefficient, which is of order erL, where a is the turbu- 
lent velocity dispersion and L is the correlation length 
of the turbulence. Observations indicate that the turbu- 
lent motions in molecular clouds are correlated on the 



scale of the entire cloud (Heyer & Brunt 2004), so L is 
roughly the size of a cloud. It we consider a species that 
is distributed with an inhomogeneity comparable to or 
greater than that of the molecular hydrogen, Hi ^ H^ 2 , 
and a composition gradient on the scale of the entire 
cloud, Hi L, then tdis ~ L/a = t CT , where t cr is the 
cloud crossing time. Thus, the time required for turbu- 
lence to smooth out a large-scale chemical gradient in a 
protocluster gas cloud is of order the cloud crossing time; 
smaller scale gradients disappear more quickly, with the 
time required to smooth out the gradient varying as the 
square of its characteristic size scale. 

Star formation in a molecular cloud also appears to 
begin in no more than a crossing time after its forma- 



tion (e.g. Tamburro et al. 2008), so there is no clear 



separation m time scales between chemical homogeniza- 
tion and the onset of star formation. This suggests 
that molecular clouds cannot be too far from homoge- 
neous when they are assembled, but that the homoge- 
nization need not be total, since small residual chemical 
gradients will be wiped out as star formation proceeds. 
Thereafter, changes in element abundances, and inho- 
mogeneities in the resulting stars, may be produced by 
supcrnovae within the cloud. In principle, even a single 
supernova is sufficient to c hange the chemi c al sign ature 



measurably. For example, De Silva et al. (2007a) find 



that scatter in Fe abundance in the HK 1614 moving 
group is roughl y 0.01 dex. The Fe conte nt of the Sun 



is ss 10 3 M Q (Anders & Grevesse 1989), so the mea- 



sured |De Silva et al. | scatter corresponds to a scatter 
in mass fraction of ~ 10 -5 . A sup ernova from a 15 Mq 
star p roduces ~ 10 _1 M Q of Fe (Woosley & Weaver 



1995[) , so if that supernova occurs inside a star-forming 
cloud that continues forming stars thereafter, the change 
in iron abundance will be measurable even if the super- 
nova ejecta are mixed with 10 4 Mq of pre-supernova gas. 
The time required for a very massive star to evolve from 
formation to explosion defines the supernova time scale, 
tsN ~ 3 Myr. Thus, while the "initial conditions" for 
cluster formation are well-mixed, we only expect star 
clusters to be homogeneous if they are assembled on time 
scales shorter than tsN- 

3.2. Chemical homogeneity 

To determine under what conditions this requirement 
is satisfied, we must compare tsN to the cluster forma- 
tion timescale tf orm . There is considerable debate about 
whether s tar cluster formation occurs in a single cross- 



Elmegreen ||2007j |Pffamm-Altcuburg & Kro upa 
or over 3 — 4 c rossing times (|Tan et al. ||2006| |rlutl &| 
Stahler |[2006] IKrumholz fe Tan 1 \'M7\ [Jeffries | liM)7| 
Nakamura fe Li ||2007[ |Matzner | | 2007p . Since the latter 
is more restrictive, we adopt it. This leaves the problem 
of estimating the crossing time. To do so, we use a con- 
venient form that expresse s it in terms of a cloud's mass 
M and column density E ( |Tan et al. ]|2006[ ) 



ing time (Elmegreen 2000| |Hartmann et al. 2001 



tc 



0.95 



M 
S3 



1/4 
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where a v i r is the cloud's virial ratio ( Bertoldi & McKee 



19921, essentially its ratio of kinetic to gravitational po 



tential energy up to geometric facto rs. Observed star- 
formi ng clouds all have a v i r ~ 1 — 2 (McKee & Ostriker 
20071, so we take a v ; r = 1.5 as typical, but £ is somewhat 



more difficult to estimate. Partly-embedded clusters de- 
tected in nearby star-forming regions using 2MASS and 
Spitzer/IRAC have surface densities of a f ew 0.01 g cm -2 



( [Adams et al.||2006 

these samples are se 



Allen et al. 2007), but because 
ected in the near-infrared they are 



strongly biased against more deeply-embedded and thus 
presumably younger clusters. Indeed, these samples are 
all dominated by stellar rather than gas mass, and thus 
have likely undergone significant expansion. Observa- 
tions of gas-dominated cluster-formi ng regions show col- 



umn densities of ~ 0.1 — 3 g cm 2 (Plume et al. 1997 



Mueller et al~| 2002; S hirley et al~l2TOl l^aundez " e j a l 



20041 IFontani et al.J 2005D; the data are summarized in 



Figure 1 of Fall et al. 12UT0 ). We therefore adopt £ = 0.3 



g cm~ 2 as a fiducial value. This is toward the low end of 
the observed range, giving a high (and thus conservative) 
estimate of the crossing time. Together with an assumed 
formation time tf OTm — 4t CI gives 



3.0 



e 

02 



-V* / M* 



V 10 4 M e 



1/4 



Myr, 



(3) 



where e = M*/M is the star formation efficiency, i.e. the 
fraction of the initial cloud mass that is converted into 
stars. Observations suggest this is typically 10 — 30% 
( |Lada k, Lada|[2003| . 

Since iform is the time required to form all the stars, 
we conclude that for clusters 10 4 Mq or smaller, almost 
all of the stars will form before the first supernova occurs 
within the cluster, so the gas will be chemically homo- 
geneous. This is probably somewhat conservative, since 
only the most massive stars go supernova in 3 Myr, and 
such a massive star is unlikely to be found in a cluster 
containing only 10 4 Mq of stars. Moreover, isN is the 
time delay between formation of a massive star and the 
resulting supernova, and a massive star is not necessar- 
ily the first to form in a given cluster. If we assume 
that massive star formation occurs on average halfway 
through the cluster formation process, then we only re- 
quire tform < 2igN to ensure there is no self-pollution, 
which corresponds to a cluster mass M* = 1.6 x 10 5 Mq. 
We therefore conclude^ that clusters smaller than 10 4 
Mq should essentially all be chemically homogeneous, 
while a significant fraction of clusters up to ~ 10 5 Mq 
will also be homogeneous. This is to be expected if the 
surface density of their more massive cluster-forming re- 
gions are higher than for the lower-mass clusters. 

3.3. Unique chemical signatures 

The analysis in §3.2 allows us to derive an order-of- 
magnitude estimate of the number of unique chemical 
signatures that are needed to identify stars with their 

6 It is interesting to note what happens in the case of globular 
clusters. These form at a density of S ~ 3 g cm -2 which leads to a 
sixfold drop in the dynamical time and a uniformity mass limit of 
> 10 7 Mq. Apart from a few light elements, globular clusters are 
chemically uniform, with few exceptions (e.g. uCen; Gratton et al 
2004). 



parent clusters. Within the Solar toroid, the total num- 
ber of tags required comes from integrating equation [T] 
over the mass limits, viz. 



dAf 
~dt 



In M. 



In M, 



d 2 N 

dM dt 



d(lnM») 



6 x 10 5 Gyr" 1 (4) 



Thus, assuming that we could select stars with an 
age similar to the Sun (t Q « 4.57 ± 0.5 Gyr), we 
would need of order 6 x 10 5 unique chemical signatures. 
Bland-Hawthorn & Freeman (2004) suggested going after 
dozens of independent elements that could be measured 
accurately enough to distinguish a "weak" line strength 
from a "strong" one. We can express this as iV c = N™ 
where Nl is the number of distinct line strengths and 
m is the number of independent elements or element 
groups. After detailed consideration of what is possi- 
ble in the context of a fibre-fed high-resolution spectro- 
graph, the HERMES project has settled on 8 indepen- 
dent groups with 5 measurably distinct abundance levels, 
i.e. N c ~ 5 8 ~ 4 x 10 5 unique chemical signatures. This 
is the target specification for the million-star Galactic 
Archaeology survey at the Anglo- Australian Telescope 
(AAT) due to commence in 2012 (Freeman et al 2010). 
The independent groups (e.g. alpha elements) include 
one or more elements that show some variation with re- 
spect to [Fe/H]. 

The volume of the potential chemical space (C-space; 
Bland-Hawthorn & Freeman 2004) accessible to the 
HERMES survey is barely adequate for our proposed 
study. But further constraints will come from stellar ages 
at least for some stars. An important population is the 
subgiants for which ages can be derived from differen- 
tial spectroscopy. In particular, a differential precision 
in surface gravity of A log g < 0.1 can be used to deter- 
mine ages to ~ 1 Gyr from isochone fitting in the sur- 
face gravity vs. effective temperature plane. There are 
some hidden assumptions in arriving at equation [4] but 
we defer a discussion of these issues to a companion pa- 
per. However the analysis does emphasize the potentially 
vast amount of information that awaits us in systematic 
high-resolution surveys of millions of stars. 



4. THE SOLAR FAMILY 

Since the Sun is the star whose chemical abundances 
we know best, it is of particular interest for the purpose 
of chemical tagging. We would therefore like to estimate 
how many chemically-identical siblings the Sun might 
have. As we have shown, this largely depends on the 
mass of the progenitor cluster. 

Attempts have been made to constrain the size of the 
Sun's parent cluster using both both dynamical and nu- 
cleosynthe tic evidence. Dynamically, Adams & Laugh- 
(2001) argue that the nearly circular orbits of the 



lin 



outer planets suggest that the Sun was born in a cluster 

with no more than M® tC \ ~ 10 4 Mq. In a larger cluster 
the Sun would likely have experienced a close encounter 
that would have left Neptune with an elliptical orbit, and 
possibly stripped it en tirely. However, as summarized by 



Megeath et al. ( 2008 ) , t his estimate is pr obably a bit low 
due to effects that [Adams & Laughlin did not include. 
Neptune probably migrated to its present position from 
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closer to the Sun over a time scale of ~ 10 Myr, so dur- 
ing this period, when the encounter rate with other stars 
would be h ighest, the Solar Sy stem presented a smaller 
target than |Adams fc Laughlin | assumed. Moreover, even 
if the orbit of Neptune had been perturbed, its eccentric- 
ity might have damped due to interactions with a plan- 
etesimal disk. These effects imply that the dynamical 
constraint is weaker than Adams & Laughlin's estimate 
Mq c i ^ 10 4 -Mq, although by how much remains quite 
uncertain. We should add that on top of these plane- 
tary d ynamics effects, there is considerable uncertainty 
in the Adams fc Laughlin calculation arising from the 
assumptions they are forced to make about the primor- 
dial binary fraction, the initial stellar density of the Sun's 
parent cluster, and how gravitationally unbound it was, 
all of which are quite uncertain. 

The nucleosynthetic constraint comes from the pres- 
ence of the daughter produc ts of short-lived radio nuclides 
in Solar System meteorites (Wadhwa et al. 2007). These 
isotopes have half-lives in the range U.I — 6.7 Myr (for 
41 Ca and 53 Mn, respectively), so their equilibrium abun- 
dances in the interstellar medium are very small, well 
below the levels inferred from meteorites. Thus a local 
and non-equilibrium origin is required to explain their 
presence. Some of the isotopes, such as 10 Be, were likely 
produce d in situ in the ea r ly Solar System b y x-ray irra- 



diation (Shu et al. 1997 Lee et al. 1998), but others 



must be the result of stellar nucleosynthesis and injec 
tion of material into the proto-Solar System, either be- 
cause that e lement is not pro duced efficiently by irradi- 
ation ( 60 Fe; Lee et al~^|l998 ), because the x-ray flux of 
young stars is insufficient to produce enough of it in sit u 



( 7 Be, 26 Al, 26 C1, and 53 Mn; |Duprat fc Tatischeff~||2008| ), 



or because timing constraints show that the element m 
question must have been injected in a single b urst rather 



than produced over an extended period ( 26 A1; Thrane et 



al. 2006 Connelly et al. 20081. These elements must 



have an extrasolar origin, which is usually taken to be 
nucleosynthesis in a nearby massive star (although see 
Gounelle &: Meibom |2008| who argue for accretion of 
DU Fe and possibly ZD A1 from the ISM). 

The need for production in a nearby massive star sets 
constraints on the size of the Sun's parent cluster and its 
chemical homogeneity, albeit quite uncertain ones. There 
are three main scenarios for how supernova-produced el- 
ements might be incorporated into the gas out of which 
meteorites containing short-lived radionuclides eventu- 
ally formed. First, a supernova explosion might have 
triggered collapse of a pre-stellar gas core a few pc away, 
leading to formation of the Sun in a gas clo ud that has 



ust b een enriched by supernova ejecta (e.g. |Boss et al 
2008). If this scenario is correct, then we have little 



rope of finding the Sun's siblings, because the supernova 
that injected the nuclides would also lead to significant 
chemical inhomogencity in the Sun's parent cluster. Ex- 
plaining the observed abundances of short-lived radio- 
nuclides requires that the pre- Solar core capture ~ I0~ 4 



of the total supernova ejecta (Meyer & Clayton 2000) 
and, as noted above in § [3j capturing even one part m 
10 5 of the supernova ejecta would be enough to change 
the iron abundance by a measurable amount. Thus stars 
that formed prior to the supernova would be chemically 
different than those that formed subsequently, and pre- 



stellar cores at different distances from the supernova 
blast wave would be enriched by different amounts. Thus 
M G . c i would be large (probably I0 4 - 10 5 Mq or more 
in order to ensure a long enough cluster formation time 
to allow a supernova), but most of the stars would be 
chemically different from the Sun. 

A second scenario is that the short-lived radionuclides 
were delivered via the winds of a Wolf-Rayet (WR) star 
rather than by a supernova ( Arnould et al.|2~006 Gaidos 
et al.|[2009 ). The WR star would not have been part of 



the Sun's cluster, but would instead have been part of an 
earlier generation of stars born in the same giant molec- 
ular cloud (GMC), and would have contaminated that 
CMC. This scenario still requires supernovae to explain 
the presence of 60 Fe, since this isotope is not present in 
significant quantities in Wolf-Rayet winds, but its rela- 
tively long half-life of 1.5 Myr, compared to 0.73 Myr 
for 26 Al and 0.1 Myr for 41 Ca, also allows the stars re- 
sponsible for producing the contamination to have been 
part of an earlier generation rather than part of the Sun's 
parent cluster. Under this scenario we therefore have no 
constraint on the size of the cluster in which the Sun was 
born, since none of the contaminating stars were part of 
the Sun's cluster. Since the contamination in this sce- 
nario happens prior to formation of the Sun's cluster, 
there should be time for mixing to homogenize the con- 
taminants and the uniformity of the cluster should not 
be compromised. 

A third scenario incorporates the injection of short- 
lived radionuclides into the Sun's protoplanetary disk 
after star formation in the cluster is complete and the 
Sun's parent core has been fully accreted, but before both 
the cluster and the Sun's protoplanetary disk have dis- 
persed (Ouellette et al. 2007). Since in this case the 
supernova ejecta mix with a ~ 0.01 Mq protoplanetary 
nebula rather than the entire 1 Mq pre-stellar core, the 
proto-Solar System need capture only ~ 10 -6 of the su- 
pernova ejecta to explain the observed abundance, and 
the chemical change induced by capturing such a small 
fraction of a supernova's ejecta is generally below the 
observationally-determined level of chemical homogene- 
ity in star clusters. Moreover, in this scenario the Sun 
should be enriched much less than the proto-planetary 
nebula by the supernova, since its small cross section 
ensures that it captures very little supernova material 
directly, and only a small fraction of the enriched gas in 
the protoplanetary disk will accrete onto the Sun rather 
than being dispersed. Thus, this scenario is still compati- 
ble with the Sun being part of a chemically-homogeneous 
cluster. 

We can estimate how large the Sun's parent cluster is 
likely to have been under the disk enrichment scenario 
based on the likelihood of a supernova occurring close 
enough to provide the observed level of enrichment and 
early enough to precede the disk's dispersal. Small clus- 
ters are unlikely sites for the formation of the Sun be- 
cause they do not produce stars massive enough to go 
supernova before the protocluster disperses, while clus- 
ters that are too large are unlikely because they have 
short re laxation times that cause the m to spread apart 
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to the degree required by observations is most likely in 
a cluster of ~ 10 4 stars, with a broad peak around that. 
(The absolute probability depends on the assumptions 
one makes about the initial cluster surface density, the 
duration of star formation, and the survival fraction of 
protoplanetary disks that come close to ionizing stars for 
at least parts of their orbits, but these do not change the 
most likely cluster size.) Given the uncertainties, we can 
state that under this scenario we expect the Sun to have 
been born in a cluster of 10 3 — 10 5 Mq in mass. 



5. IN SEARCH OF THE SOLAR FAMILY 

In this section, we discuss basic parameters of the Solar 
family (e.g. its orbit properties) and of the Galaxy disk 
to allow us to simulate the projected distribution of this 
particular group. From these simulations, we are led to 
possible survey strategies for finding group members. 

5.1. Solar parameters 

We know more about the Sun than any other star 
which makes it the best template for setting up a chemi- 
cal tagging experiment. The Solar atlas of spectral lines 
includes 25,000 lines covering 66 elements (Moore, Min- 
naert & Houtgast 1966; Wallace, Hinkle & Livingston 
1998). We know a great deal about the Sun, including its 
absolute abundances in all elements, its age, mass, lumi- 
nosity, and its surface temperature. We know the depth 
of the convection zone, its central temperature, internal 
structure, and the strength of its magnetic field. We 
have some understanding of its surface activity, surface 
structure and complex rotation. In principle, with access 
to a vast database of echelle spectroscopy, it should be 
possible to identify the Solar family more easily than any 
other disrupted star group, at least for Solar type stars. 
This task will be greatly aided in an era of ESA Gaia 
when candidacy for the Solar family can be restricted by 
complete phase space information (see §7). 

We assume the Sun was born in a dense cloud which 
gave rise to 10 4 stars to within some factor. A star clus- 
ter of this mass may have diffused within the last few 
billion years, but initially we shall assume that the star 
cloud has dispersed uniformly within a torus along the 
Solar circle. Since proof of membership is obtained from 
chemical tagging (via heavy element abundances mea- 
sured from high-resolution spectroscopy), we examine 
how many stars exceed apparent magnitude limits ap- 
propriate to modern telescopes, i.e. apparent magnitude 
limits in the range B = 14 to B = 20. These are the likely 
flux limits for a wide-field, multi-object spectrograph on 
any existing or foreseeable telescope. 

We compute the present-day colour-magnitude dia- 
gram (CMD) for the putative Solar family using the pro- 
gram StarFIStQ We specify BASTI isochrones with the 
revised Solar mctallicity Z = 0.0198 (corresponding to 
[Fe/H]=0.06) for the known age of the Sun, 4.57 Gyr. A 
Salpcter IMF (logarithmic slope —1.35) is assumed with 
a 25% binary star fraction. Our total population is 10 4 
stars with a minimum stellar mass of 0.1 Mq and max- 
imum mass of 100M Q . StarFISH returns the data in 
pixelated form, for which the colour range sampling is 

7 See marvin.as.arizona.edu/~jharris/SFH 
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Fig. 1. — Predicted B-V and B-R colour-magnitude diagrams for 
the Solar family; the parallel main sequence arises from unresolved 
binary stars. The total population of simulated stars is 10 4 ; the 
greyscale indicates the relative number of stars in bins of 0.05 along 
both axes. 



(B-R and B-V): to 3.0 ( in steps of 0.05); and magni- 
tude range (R and V): 10 to -3.0 (in steps of 0.05). The 
resulting CMDs are shown in Fig. [I] 

We now consider the Sun's position in a cylindrical 
frame. The radius of the Solar circle is close to 8.0 kpc 
(Reid 1999). Gies & Helsel (2005) have determined a use- 
ful orbit for the Sun within Dehnen & Binney's (1998) re- 
alization of the Galactic potential. The radial excursion 
keeps the Sun confined to the radial range (R : 7.8 — 8.6) 
kpc essentially consistent with an epicyclic approxima- 
tion for a nearly circular orbit in an axisymmetric grav- 
itational field (Lindblad 1959). If <jr is the velocity dis- 
persion (in km s _1 ) in the radial direction, and r\ is the 
rms radial amplitude of the epicyclic motion in kpc, then 
<7r = 36ry. For the old stellar disk, an is about 40 km 
s" 1 (Nordstrom et al 2004). 

Thus the Sun is presently near pericentre and its guid- 
ing centre is at 200 pc beyond the Solar circle. Gies 
& Helsel determine an excursion period close to the or- 
bit period and a frequency that is essentially identical 
to the value deduced from Oort's constants (36.7 ± 2.4 
km s _1 kpc -1 ). While the Sun was long thought to be 
substantially offset from the disk plane, the Hipparcos 
astrometric satellite has convincingly demonstrated that 
its present position is 8±4 pc (Holmberg, Flynn & Lindc- 
gren 1997). Gies & Helsel determine a vertical oscillation 
amplitude of 140 pc roughly a factor of three less than 
what is expected for the old stellar disk (§5.3). Thus 
for the Solar toroid, we adopt a rectangular cross-section 
with radius 8.2 kpc, a radial thickness 800 pc and a ver- 
tical thickness of 140 pc (scale height 70 pc). In §6, we 
refer to this as Model 1. 

5.2. Galactic extinction 

In any projection of a dispersed stellar group, we must 
include the effects of interstellar extinction along the line 
of sight. We investigate this problem in two ways which 
gave similar results. First, we generate a fractal medium 
on the scale of the Galactic disk. This is assumed to be 
a flat disk with a vertical thickness of 500 pc and with 
a total radial extent of 20 kpc. After Fischera & Do- 
pita (2005), we adopt a lognormal dust-gas distribution 
consistent with large-scale turbulence in a compressible 
medium. The probability density distribution P of the 
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local density p in the medium is given by 



P(p) = -P(lnp) = 



1 

/ Z7T(TpP 



,-*V(2<t 2 ) 



(5) 



where x = In p — In p a , cr 2 , is the variance and p is the 
mean of the lognormal density distribution. If dust is a 
passive contaminant in this medium, then its local den- 
sity distribution is the same as the gas density (Warhaft 
2000). There are several methods for implementing such 
a medium: one of the most approachable is given by 
Elmegreen (1997) but a more rigorous treatment uses the 
Fourier transform domain (Sutherland & Bicknell 2007). 
Fischera & Dopita (2005) show how to normalize the syn- 
thetic medium to account for local extinction by the ISM 
at optical wavelengths. 

Our second method is to utilize the IRAS DIRBE ex- 
tinction maps derived by Schlegcl, Finkbcincr & Davis 
(1998F] At the present time, there is no 3D model of 
how the dust is distributed along each and every line of 
sight. We therefore adopt a statistical approach. The 
dust is assumed to be distributed uniformly along each 
line of sight. The total extinction E(B-V) is taken di- 
rectly from DIRBE maps with a conversion to A# taken 
from Table 2 in Schlegel et al (1998). The columnar ex- 
tent of the dust distribution is assumed to be identical 
to the disk used in the above method. Our approach en- 
sures that if we were to measure the extinction along each 
line of sight with an essentially infinite number of model 
stars, our synthetic extinction map would still match the 
complexity of the observed extinction map. 

5.3. Radial and vertical diffusion 

When a star experiences deflections along its orbit, the 
phase space density around it changes with time. It is 
well known that older stars have larger velocity disper- 
sions than younger stars and, indeed, the study of stellar 
diffusion has a long history (e.g. Spitzer & Schwarzschild 
1951; 1953). Binney & Tremaine (2008; hereafter BT) 
give the Schwarzschild equation which describes the ve- 
locity distribution of stars in the solar neighbourhood. 
At a fixed age, the stars are described by a velocity ellip- 
soid with velocity dispersions o\r, and er 2 with respect 
to the radial, azimuthal and vertical axis respectively. 
The axis ratios are roughly constant as a function of stel- 
lar age but the volume of the ellipsoid grows with time, 
at least for a few Gyr (Soubiran et al 2008; cf. Aumer & 
Binney 2009). 

A crude approach to understanding this behaviour is 
to assume kinematic diffusion using coefficients that are 
constant and isotropic (Wielen 1977; cf. Rosenbluth et 
al 1957). This model, which is straightforward to im- 
plement within a numerical code, considers small veloc- 
ity perturbations Avi and radial perturbations Ar^ that 
have a mean of zero when averaged over all perturba- 
tions (denoted by the subscript i) over some period of 
time. But the rms value need not vanish such that it 
increases with time At, 



\2\ 



= S v At 



(6) 



8 These authors supply a code for accessing t he DIRBE maps 
at http://www.astro.princeton.edu/~schlcgcl/dust which we have 
adapted tor compilation under MAC USX: a copy can be obtained 
from jbh@physics.usyd.cdu.au. 



where 8y is the velocity diffusion coefficient (Sy = 
8v(8r^4>,8z))- This diffusion process also acts on the 
stellar positions such that 



(£(Ar 4 ) 2 )H]T(A^-^)) 2 



(7) 



where the perturbation at time £j led to the new posi- 
tion at time t. For the observed vertical oscillation and 
in-plane epicyclic frequencies, ui z and ujr^, the rms ex- 
cursions are then obtained in Galactocentric radius R 
and vertical height z from 

((AR) 2 )= f S R sinc 2 u Rlj> {t-t) 

J O 

+ §4,(1 - cosuj R ^t - t')) 2 /(2B Q ) 2 dt (8) 
((Az) 2 ) = / S„ smc 2 uj z {t-t)dt (9) 



for whic!^]-B© is Oort's constant. Wielen (1977) con- 
sidered three empirical relations for <5y: (i) a constant 
coefficient; (ii) a velocity dependent coefficient; (ii) a 
velocity-time dependent coefficient. We solve for each 
of these and present the results in Fig. 2. The vertical 
dotted line indicates the age of the Sun and, as shown 
here, the characteristic scale height for a middle aged 
star is 200-250 pc. The rms radial spread is ~2 kpc, i.e. 
about a factor of three larger than realized by Gies & 
Helsel (2005). 

But the velocity of most stars changes significantly 
along an orbit which argues for a more sophisticated 
treatment. Indeed, Binney & Lacey (1988) considered 
how stars diffuse in integral space and showed that the 
diffusion coefficients vary from point to point. They con- 
sidered Spitzer & Schwarzschild's proposal of random 
scatterings by giant clouds and found that, with the as- 
sumption of energy equipartition, o~ z should lie between 
(70 and o~ R . But this is not what is observed. The ratio 
of the last two quantities is fixed by the epicyclic and 
azimuthal frequencies in the solar neighbourhood (BT; 
equation (8.117)), but a z is observed to be smaller than 
both a R and cr^ (Nordstrom et al 2004). Ida, Kokubo 
& Makino (1993) and Shiidsuka & Ida (1999) challenged 
Binney & Lacey's equipartition argument in a revised 
analysis and find that a z is the smallest component of 
the dispersion, as observed. This behaviour has been 
verified by Sellwood (2008) who demonstrates that the 
anisotropic distribution of the (more distant) massive 
clouds determines the shape of the velocity ellipsoidj^] 

Some early papers (e.g. Carlberg 1987; Jenkins & Bin- 
ney 1990) appeal to transient spiral arms to assist with 
the vertical heating. But because the vertical oscillation 
frequency is higher than the radial epicyclic frequency, 
spiral scattering couples to vertical scattering only very 
weakly. Thus, the processes that drive vertical scattering 
arc distinct from churning. In a recent twist, Schonrich 
& Binney (2009) investigate the idea that at least some 
of the vertical dispersion is assisted by the churning pro- 
cess; here, stars that have migrated from the inner disk 
retain their higher vertical dispersions when migrating 

9 The sine cardinal function is defined as sinc(ir) = x~ 1 shirr. 

10 Presumably this is strong evidence against the recurrent claim 
of dark matter halos consisting of massive black holes. 
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outwards. But whether non-circular orbits can be effi- 
ciently transported in this way remains unclear. 

Aumer and Binney (2009) define an empirical velocity 
dispersion a 2 for a known stellar distribution TV of a given 

age t as 

^ = a 2 (t)(dN/dt) dt 
J*°(dN/dt) dt 
where the vertical component is given by 

-* (t) = 2 °(^r- (u » 

We use updated values for the coefficients (Binney 2010) 
where the age of the disk is to = 10 Gyr and for which 
ti = 0.1 Gyr. For our application, t = t a is the age (4.57 
Gyr) of the Solar Family. 

The vertical dispersion can be related to a physical 
scale, i.e. a characteristic scale height or disk thickness, 
if the vertical component of the potential & z is known 
since \a 2 z ~ |$ z | (Kuijken & Gilmore 1991). A more 
rigorous treatment includes information from the vertical 
distribution function f z , but the existing data do not 
provide a consistent story (see Binney 2010). Suffice it 
to say, the mean scale height for solar-type stars is close 
to z a s=y 250 pc (Wielen 1977) which is confirmed by the 
more sophisticated analysis. 

5.4. Radial transport 

In the Widen picture, star clusters dissolve and grad- 
ually diffuse into the background galaxy. Stellar popula- 
tions as an ensemble experience heating. The Scllwood- 
Binney mechanism is quite different in character where 
perturbations allow the star to migrate in integral space 
without heating the disk. During interaction with a sin- 
gle steady spiral event of pattern speed flp, a star's en- 
ergy (E) and angular momentum (J) change while it 
conserves its Jacobi integral. In the (E, J) plane, stars 
move along lines of constant ij = E — Op J. For tran- 
sient spiral perturbations, the star undertakes a random 
walk in the (E, J) plane, deflected by a series of presum- 
ably uncorrelated spiral arm events which occur on time 
intervals of 500 Myr or so. (The impulse acting on indi- 
vidual stars lasts for a fraction of this interval.) So the 
Solar family would have experienced about ~10 of these 
events in the vicinity of co-rotation if only a fraction of 
the stars are moved in each event. 

Substantial variations in the angular momentum of a 
star are possible over its lifetime. A single spiral wave 
near co-rotation can perturb the angular momentum of 
a star over a broad distribution in J z with a tail up to 
50%. The star is simply moved from one circular orbit 
to another, inwards or outwards, by up to 2 kpc or so. 
Conceivably, the same holds true for a star cluster except 
that most are dissolved in much less than an orbit period. 
Given that there exists a long tail of very old star clusters, 
it is interesting to note that some of these could have 
been picked up and moved to safety away from the major 
disturbances. 

Support for the Sellwood-Binncy mechanism (churn- 
ing) comes from recent N-body simulations that allow 
for a steady accretion of gas from the outer halo (Roskar 
et al 2008). These relatively "cold" simulations reveal 



some of the secular processes that may be acting over the 
lifetime of the disk. Averaged over the entire disk, these 
authors find that the rms radial excursion at t = 10 Gyr 
is ((AR) 2 )"- 5 w 2.4 kpc and (AR) w kpc. For the out- 
ermost disk stars, the simulations show that (AR) rts 3.7 
kpc, several times larger than their epicyclic radial ex- 
cursion and our derived value of r] (§5.1). 

In order to understand the possible influence of churn- 
ing on the Solar family, we envisage that at least one 
event affected the Solar family, i.e. those that had the 
right pattern speed to move the members from wherever 
they were born to where they are today. These events 
need not have treated all family members the same way. 
In all likelihood, the Solar family was already dispersed 
azimuthally at the time of the first event. Depending 
on their orbital phase at the time, relative to the spiral 
event, some may have gone in and others out, and others 
may not have been much affected. There may have been 
subsequent events with the right pattern speed to move 
some members of the Solar family from their new radii 
to somewhere else. 

In order to model the impact of churning on the Solar 
Family, we use the illustrative formula provided by SB02. 
The probability that a star formed at radius R migrates 
to R at time t is given by 

P(R \R) = (2^)-°^ (- {R °~* ? ) (12) 
for which 

<7* =iZ o ((0.16) 2 + (0.4) 2 -^) - 5 (13) 

to 

Using this simple formalism, we determine that the Solar 
Family could be dispersed over an annulus ~ 5.5 kpc. 

Churning is a symmetric process enforced by angular 
momentum conservation. But a family of stars in the so- 
lar neighbourhood can all move inwards together, or out- 
wards, or even asymmetrically depending on how much 
diffusion a star cluster has experienced and the conse- 
quent phase lags with respect to the disturbing force 
(e.g. SB02, Figs. 4 and 13). Asymmetries can take 
place near inner (ILR) or outer Lindblad (OLR) reso- 
nances in the sense that ILRs scatter inwards and OLRs 
scatter outwards. Interestingly, both phenomena appear 
to have taken place in the solar neighbourhood in re- 
cent times. Dehnen (2000) has invoked an outer Lind- 
blad resonance driven by the central bar to explain the 
Hercules Stream. Sellwood (2010) finds good evidence 
from the Hyades Stream for an inner Lindblad resonance 
which presumably arose from a recent spiral disturbance 
that extended far out into the disk. This is evidence for 
the phenomenon of transient spiral waves that propagate 
throughout the Galaxy due to large-scale disturbances to 
the disk. In the simulations of Roskar et al (2008), there 
is clear evidence of asymmetry of outward vs. inward 
transport at a radius of 8 kpc. But this reflects the lim- 
ited extent of the star formation beyond this radius. 

To investigate the signature of a skewed distribution, 
we derive an empirical function that has a skewed Gaus- 
sian cross-section in the radial direction. We assume a 
coeval population that is distributed around the Solar 
toroid The Gaussian function has the form 

G(R-R )=p exp[-ln2-g(R-R )} (14) 
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Fig. 2. — Stellar diffusion in radius (solid curves) and vertical 
to the disk (dashed curves) as a function of cosmic time. The 
different time-dependent relations (see text) were obtained from 
three different diffusion coefficients first derived by Wielen (f977). 
The vertical dotted line indicates the age of the Sun. 
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Fig. 3. — The skewed gaussian function for four values of skew- 
ness: s = (solid), s = 0.3 (dotted), s = 0.65 (dashed), s = 1.0 
(dot-dashed). These plots are normalized to the same peak den- 
sity but can be normalized to the total area under the curve (see 
text). This function is used to describe the cross-section of the 
Solar toroid for different amounts of radial diffusion. 

where g(x) — (2x/J r g ) 2 . The quantity T g is the 
FWHM (i.e. natural width) of the Gaussian toroid. 
The cross-sectional area of the Gaussian toroid, A g — 
(J r g /2)(ir/ In 2) ' 5 , allows the total volume of the toroid 
to be easily derived. 

This is transformed to a skewed function by replacing 
g(x) with h{x) given by 

h{x) = (ln[l + 2sx/T h ]/s) 2 (15) 

where s is a measure of skewness and x > — J r / l /(2s). 
This has the nice property that h{x) —> g{x) in the limit 
that s — > 0. We can therefore write 

H{R-R ) =p exp[- In 2- h(R - R )} (16) 

In Model 3 below (§6), the variable s is time depen- 
dent. The quantity T h is the FWHM (i.e. width) 
of the skewed Gaussian toroid which is a function of 
s. The cross-sectional area of the skewed toroid Ah = 



(JV2)(7r/ln2) - 5 exp[s 2 /(41n2)] allows the total volume 
of the toroid to be approximated or for the total area 
to be renormalized to the gaussian cross-section with 
the same FWHM. In Fig. [3j we plot the function for 
s = 0,0.3,0.65, 1. Radial profiles with s ~ 0.5 resemble 
the slightly skewed histograms presented in Fig. 13 of 
SB02. We emphasize that while certain annuli within 
the disk may exhibit this behaviour, churning conserves 
angular momentum when summed over all sections of the 
disk. 

6. RESULTS 

6.1. Distribution in cylindrical coordinates 

In Fig. HJ we show four possible distributions for the 
Solar family today. The model parameters are inspired 
by different published studies: 

1. toroid parameters derived from orbit modelling by 
Gies & Helsel (2005) - the radial and vertical thick- 
nesses (FWHM) are 800 pc and 140 pc respectively 
(see §5.1); 

2. toroid parameters derived from equations 10-13 - 
the radial and vertical thicknesses (FWHM) are 5.5 
kpc and 500 pc respectively (see §5.3 and §5.4); 

3. toroid parameters derived from Wielen et al (1996) 
and equations 10-13 where the stars are born at a 
radius of R — 6 kpc - the radial and vertical thick- 
nesses (FWHM) are 4.1 kpc and 500 pc respectively 
(see §5.3 and §5.4); 

4. toroid parameters to illustrate a skewed distribu- 
tion - the vertical thickness (FWHM) is 500 pc 
and the radial spread is asymmetric (s = 0.5 at the 
present epoch; see §5.3 and §5.4). 

Note that in polar coordinates, the stellar density must 
be weighted by a factor of R/R to ensure a radial density 
that conforms to the distributions in Fig. [3j 

6.2. Distribution in Galactic coordinates 

In Figs. 5-8, we illustrate the distribution of Solar 
siblings in Galactic coordinates for the four models in 
Fig. 4. Only stars above a magnitude limit B = 18 are 
shown. The blue box in each figure aids a comparison 
between the distributions. This region falls inside of the 
LSR apex; we refrain from illustrating the mirror field 
that is symmetric about (£,b) = (0,0). The distribution 
in longitude I is quite different for each of the models. 

Note that the blue box in all figures is centred at I = 
75° rather than in the direction of the Sun's motion. 
This is because we are witnessing the integrated effect of 
the path length as the Solar torus bends inwards. This 
is particularly apparent for Models 2 and 4. The high 
projected density of stars in Model 3 reflects the higher 
surface density when projected onto the disk compared 
to Models 2 and 4. 

The stars in Model 1 are characteristically nearer to 
the Sun; interestingly, useful numbers of these stars are 
predicted to fall within a B — 14 survey limit. Unsurpris- 
ingly, the stars in Model 1 are also concentrated closer to 
the disk than in the other models which may make them 
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Fig. 4.— (a) Model 1: parameters from Gies & Helsel (2005); (b) 
Model 2: parameters from Sellwood & Binney (2002); (c) Model 
3: parameters from Wielen et al (1996); (d) Model 4: parameters 
to illustrate asymmetries observed in numerical simulations. The 
white circle shows the Sun's present position. 



harder to observe in practice due to crowding. The offset 
position of the Sun in relation to the thin torus (see Fig. 
4(a)) leads to the asymmetric distribution of siblings on 
the sky at both low and high galactic latitude, particu- 
larly in directions that look outside of the Solar Circle 
(£ : 90° - 270°). 

The cumulative counts are presented in Fig. [9] (all 
sky) and tabulated in Table 1. The counts specific to 
the box region {£ : 60° — 90°) and neighbouring fields 
are also given in Table 1. We refer to longitude zone 
(£ : 90° — 120°) as the outer region and longitude zone 
(£ : 30° — 60°) as the inner region. The trends observed 
in Table 1 between the inner and outer regions largely 
reflect what we see in Figs. 5-8. Models 2 and 3 have 
a substantial inner population with a small fraction of 
outer stars in the case of Model 3. Models 1 and 4 have 
a substantial outer population which undergoes a rever- 
sal at faint magnitudes. Prima facie, models 2 and 4 
look similar but the behaviour exhibited by the magni- 
tude counts in Table 1 is different. 

Our models are simply to illustrate the prospect of di- 
rectly measuring the impact of secular evolution within 
the disk. We stress here that a more sophisticated ap- 
proach to dust extinction is needed to be sure of the num- 
ber counts, particularly at the bright magnitude limit. 
This is particularly true for Model 1 where the projected 
stars lie close to the Sun. We refer the reader to a 
brief commentaries on this issue by Robin (2009) and 
by Aumer & Binney (2009). 

6.3. Survey strategy 

As we saw in the last section, the distribution in z and 
the centroid position in longitude are powerful discrim- 
inants for determining the Sun's origins and the impor- 
tance of in-plane diffusion and transport. We now briefly 
consider possible survey strategies for four different ap- 
parent magnitude limits. We restrict our discussion to 
the B band because this spectral region provides the rich- 
est information for chemical tagging (Bland-Hawthorn & 
Freeman 2004). For all of the strategies below, it should 
be possible to achieve the required number of chemical 
tags derived in §3.3. 

In Fig. [9j we show four possible B magnitude limits 
(indicatedDy the vertical lines) for high-resolution spec- 
troscopy. The canonical values are as follows: 

• B=14 is an appropriate limit on a 4m class wide- 
field spectrograph operating in 2" seeing. This is a 
magnitude brighter than the effective limit of the 
HERMES instrument (V~14) that sees first light 
at the AAT in 2012. 

• B=16 is a suitable limit for a HERMES-style in- 
strument carrying out a shallow survey on an 8m 
class telescope in 1" seeing. 

• B=18 is an appropriate limit for a HERMES-style 
instrument carrying out a deep survey on an 8m 
telescope in 1" seeing. 

• B=20 is a suitable limit for a HERMES-style in- 
strument undertaking a million-star survey on an 
extremely large telescope (ELT) in 1" seeing. 

Table 1 indicates the number of Solar Family stars ex- 
pected to be found in a million-star survey within the 
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blue box region for different magnitude limits for Models 
1-4. Notice that in each of the models, the maximum 
number of Solar siblings is detected for surveys with a 
limiting magnitude of B = 16-18. This peak is a re- 
sult of the competition between two effects. For limiting 
magnitudes of B = 16 or less, only the brightest few So- 
lar siblings are included in the sample, so increasing the 
limiting magnitude increases the number of Solar sib- 
lings found. However, the total number of stars in the 
survey field that are not Solar siblings also increases as 
one goes to fainter limiting magnitudes, and thus as the 
limiting magnitude B increases, the fraction of surveyed 
stars that are Solar siblings begins to decrease. Beyond B 
= 18, this decrease in the fraction of stars that are Solar 
siblings outweighs the increase in the absolute number of 
Solar siblings brighter than the magnitude limit, and the 
expected number of Solar siblings discovered in a survey 
of fixed size therefore decreases as the magnitude limit 
gets fainter. Thus the optimal magnitude limit of B = 16 
to B = 18 is determined by a combination of the distribu- 
tion of Solar family B magnitudes, illustrated in Fig. [9j 
and the distribution of B magnitudes for all stars in the 
Galaxy, which rises sharply above B = 18. The brighter 
limit is to be preferred because of the higher signal and 
because of the reduced contamination from fainter stars 
close to the Galactic plane. 

A reasonable strategy is to survey the blue box region 
and control fields of similar size on the other side of the 
LSR apex (see columns 2 and 4 in Table 1). The longi- 
tude test may be relatively free from identification errors 
although a possible complication is the colour selection 
of the input catalogue. The more numerous faint stars 
with lower apparent magnitudes will be more reddened 
than brighter stars on average since they will be stars 
that are embedded deeper in the torus. However, stars 
that are far from the Solar torus in longitude will be less 
reddened. A selection criterion will need to take this 
into account. In an era of ESA Gaia, with the advent 
of proper motions and 6D phase space information for 
millions of stars, this test can be conducted with higher 
sensitivity in a physical coordinate frame. 

We refrain from a detailed analysis of survey strategies 
at this time because there are many issues to consider in 
an observing campaign. The model distributions them- 
selves impose their own challenges. If the Gies & Helsel 
solution for the Solar orbit is correct, the vertical con- 
centration will keep the stars much closer to the Galactic 
plane. This has two problems: higher galactic extinction 
on average and stellar crowding. In fact, for bright stars 
this may not be a serious problem even in the plane. The 
APOGEE survey is targetting red clump stars and red 
giants down to H«13 between b = —15° and b = 15° 
for the entire longitude range observable from the north- 
ern hemisphere (Allende Prieto et al 2008). This survey 
makes use of the SDSS 2.5m telescope under natural see- 
ing conditions. The chemical information appears to be 
limited at the present time although this is currently un- 
der investigation (J. Lawler, personal communication). 
The H-band limit is equivalent to roughly B«16 and 
so somewhat deeper than the HERMES survey at the 
AAT3.9m telescope that will target spectral windows in 
the V and R bands. 

In the search for the Solar family, dense stellar fields 
can be an advantage with wide-field AO systems and 



multi-object spectrographs. In particular, one can com- 
fortably envisage high-resolution IR spectrographs oper- 
ating with multi-conjugate adaptive optics in the next 
decade, either on 8m class telescopes or on ELTs. The 
"echidna" positioning technology permits 4000 stars to 
be observed simultaneously as was originally proposed 
WFMOS instrument (Colless 2005). Another related ac- 
tuator technology provides the basis for the LAMOST 
project (Zhang & Qi 2008). These new developments 
can be extended to even larger formats on wide-field fo- 
cal planes. It is reasonable to suppose that such a facility 
linked to an AO system will become available in the next 
decade because inter alia some of the most ancient stars 
are expected to have settled to the inner bulge (Bland- 
Hawthorn & Peebles 2006; Tumlinson 2010). 

7. DISCUSSION 

Within the domain of Solar physics, planetary science 
and geophysics, the formation and evolution of the So- 
lar nebula is of paramount importance. This arena has 
been revitalised by the discovery of exoplanets around 
nearby stars. In certain important respects, a search for 
the Solar family is the "missing link" between top-down 
models of Galaxy disk formation, and bottom-up mod- 
els of Solar system formation. The search for the Sun's 
siblings is the first step to placing its progenitor cloud in 
the context of the interstellar medium and the Galaxy as 
a whole. 

We have seen that a reasonable number of the Sun's 
siblings can be identified in a targetted million-star sur- 
vey, at least in principle. Given that strong resonances 
have clearly influenced the solar neighbourhood at dif- 
ferent times (Dehnen 2000; Sellwood 2010), it would be 
surprising if the Solar family has escaped significant mi- 
gration into the surrounding disk. Even more surprising 
would be the failure to identify a single family member. 

Indeed, we would argue that chemical tagging, or a re- 
lated technique, is demanded by secular processes like the 
Sellwood-Binney mechanism. If this process is as strong 
as implied by recent simulations, then most stars may 
not have formed in situ, which is an axiom of chemical 
evolution models. In such an event, a certain degree of 
reconstruction is called for although much can be learnt 
from the higher moments of the chemical space. The fact 
that ancient star clusters like Collinder 261 (De Silva et al 
2007) have survived this process may indicate that these 
were picked up and transported early in their history out 
of harm's way. 

In §3.3, we argued that full-blown chemical tagging re- 
quires N c ~ 6 x 10 5 unique tags per Gyr which may be 
within reach of the HERMES survey in its present in- 
carnation. But here the impending impact of the Gaia 
astrometric mission becomes clear. The Sun's orbit is 
one of the more circularized of nearby dwarfs (Robles 
et al 2008). Relative to other stars of similar age, its 
peculiar velocity is in the lowest 8% of a Maxwellian dis- 
tribution, consistent with the orbital eccentricity derived 
by Robles et al (2008). The Sellwood-Binney "churning" 
process does not impart any random energy (i.e. heat- 
ing). Thus, accurate phase space information from Gaia 
once available can allow us to preselect orbit families, 
for example, stars with circularized orbits. The bright 
end of this distribution (say V<14) then provides the 
primary input catalogue for a chemical tagging survey 
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where Gaia is expected to yeld its most accurate phase 
space information. 

The use of phase space information to aid the recon- 
struction of ancient star clusters is fraught with prob- 
lems at the present time. The success or failure of this 
approach relies on the assumption that the Sun's low pe- 
culiar velocity will be shared by all the other stars in 
the Solar family. Whether this is true or not depends on 
whether the Sun has a low peculiar velocity simply be- 
cause random diffusion processes just happened to put it 
at the low velocity tail of the distribution, or because of 
some property of its parent cluster that will be shared by 
all stars that were born in it. This uncertainty is partly 
reflected in our choice of three different models for the 
Solar family distribution based on very different assump- 
tions. Clearly, some knowledge of the distribution of the 
Solar family will go a long way to resolving these issues, 
and the importance of secular processes in general. 

In our attempt to recover the Solar family, we will 
identify many thousands of distinct stellar groupings in 
chemical space. In our companion paper, we derive the 
expected number of distinct stellar groupings and the 
inferred progenitor mass distribution (Sharma et al, in 
preparation). Of these long-dissolved star clusters, some 
will have much better representation, and many will 
have too few members for a reliable identification. We 
can construct the colour-magnitude diagram of candi- 
date groupings to see if they are consistent with a coeval 
population. 

Just how the reconstructed ancient cluster mass distri- 
bution will help to unravel the processes of Galaxy forma- 
tion will be investigated in future papers. In large part, 
it will depend on the detailed structure of the C— space, 
of which very little is known. There is widespread evi- 
dence of anomalous chemical signatures in many discrete 
stellar systems (De Silva et al 2008). Just how these sig- 
natures relate to the multi-dimensional topology of the 
chemical space will require much larger high-resolution 
spectral surveys to answer properly. 

At the present time, it is unclear what constitutes 
the true building blocks of galaxies, whether these are 
early-stage dwarf galaxies (i.e. the progenitors of those 
that orbit the Galaxy) or simply gas clouds confined by 
infalling dark-matter halos. But the newly discovered 
"ultra-faint" dwarf galaxies (Simon & Geha 2007) may 
also be an important constituent of the formation pro- 
cess. These may carry unique chemical signatures of the 
first stars. The faint end of the dwarf galaxy luminosity 
function can only be studied within a few megaparsecs of 
the Galaxy. Many of these systems may have already dis- 
solved within the Galaxy to be revealed in future "chem- 
ical tagging" and astrometric surveys. 

It is clear that we are only now beginning to under- 
stand just how much information is locked up in stellar 
atmospheres. In the ESA Gaia era, we will presumably 
identify tens of thousands of stellar groupings in phase 
space or in integral space. Many of these may be spurious 
dynamical groupings created by spiral arm density waves 
(De Simone, Wu & Tremaine 2004; Quillen & Minchev 
2005). But many of these may be disrupted ancient star 
clusters, some formed within the Galaxy, some external 
to it. Some of these systems may even carry the hall- 
marks of the first generations of stars in the early uni- 
verse (Tumlinson 2010). Once again, chemical tagging 



90° 




-90" 



Fig. 5. — Model 1 distribution of Solar family stars brighter 
than B=18 (see Fig. 4(a)) in Galactic coordinates as seen from 
the Sun. Note that the stars are distributed symmetrically about 
the LSR apex and antapex in longitude, quite unlike what is seen 
in the other models. Of the four models presented in Figs. 5-8, 
this model is the most susceptible to how we treat the local dust 
distribution (cf. Robin 2009). The blue box is the survey region 
discussed in the text. 
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Fig. 6. — Model 2 distribution of Solar family stars (see Fig. 
4(b)) in Galactic coordinates brighter than B=18. The blue box is 
the survey region discussed in the text. 
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Fig. 7. — Model 3 distribution of Solar family stars (see Fig. 
4(c)) in Galactic coordinates brighter than B=18. Most of the 
stars fall "inside" of the galactic latitude (projected radius) of the 
blue box, discussed in the text, with few stars outside. 

will provide key information on the likely origin of these 
low-mass structures. 

If the chemical tagging technique is ultimately success- 
ful, it will provide rigid constraints on the importance of 
the Sellwood-Binney mechanism and, by implication, on 
the importance of secular processes and the lifetime of 
spiral arms. These are fundamental properties of galax- 
ies that must be understood if we are to have any hope 
of claiming a detailed physical understanding of galaxy 
evolution. 
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Fig. 8. — Model 4 distribution of Solar family stars (see Fig. 
4(d)) in Galactic coordinates brighter than B=18. Most of the 
stars appear highly localised in two clouds each covering about 
300 sq. deg. as illustrated by the blue box at £ = 60° to i = 90°. 
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apparent B magnitude 

Fig. 9. — Stellar number counts for (a) Model 1 (b) Model 2 (c) 
Model 3 (d) Model 4; see Fig. 4. The black histogram shows the 
number of stars for each bin in apparent B magnitude before dust 
correction, red shows the number of stars in each B magnitude bin 
after dust correction, and blue shows the cumulative distribution of 
stars (corrected for extinction) with apparent magnitude brighter 
than B mag. The vertical dashed lines indicate the magnitude 
limits for four different surveys — see text. 
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TABLE 1 

Cumulative B magnitude-limited star counts for the 

MODELS ILLUSTRATED IN FlG. 4. COLUMN 1 IS THE APPARENT B 
MAGNITUDE LIMIT. COLUMN 2 GIVES THE STAR COUNTS INSIDE OF 
THE BLUE BOX CENTRED AT (£,b) = (105°, 0) FOR MAGNITUDES 
BRIGHTER THAN THE B LIMIT IN FlGS. 5-8. COLUMNS 3 AND 4 
ARE THE SAME AS COLUMN 2 BUT CENTRED AT (I, b) = (75° , 0) 

and (£,b) = (45°, 0). Column 5 gives the star counts for the 

ENTIRE SKY BRIGHTER THAN THE B LIMIT. COLUMN 6 IS THE 
PREDICTED NUMBER OF SOLAR SIBLINGS IN A MILLION STAR 
SURVEY DOWN TO THE B LIMIT. 
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